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, ^ ' Abstract. Clear signs of excess dileptons above the known sources were found at the 

SPS since long. However, a real clarification of these observations was only recently 
I achieved by NA60, measuring dimuons with unprecedented precision in 158 A GeV 

■ In-In collisions. The excess mass spectrum in the region M<1 GeV is consistent 

with a dominant contribution from tt+tt^ ^ p — > fi^ annihilation. The associated 
^ ' p spectral function shows a strong broadening, but essentially no shift in mass. In 

^ . the region M>1 GeV, the excess is found to be prompt, not due to enhanced charm 

^ I production. The inverse slope parameter Tcff associated with the transverse momentum 

^ ' spectra rises with mass up to the p, followed by a sudden decline above. While the 

^ . initial rise, coupled to a hierarchy in hadron freeze-out, points to radial flow of a 

hadronic decay source, the decline above signals a transition to a low-flow source, 
presumably of partonic origin. The mass spectra show at low transverse momenta the 
^ I steep rise towards low masses characteristic for Planck-like radiation. The polarization 

of the excess referred to the Collins Soper frame is found to be isotropic. All 
observations are consistent with the interpretation of the excess as thermal radiation. 
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Q ■ Dileptons are particularly attractive to study the hot and dense QCD matter 



O 



formed in high-energy nuclear collisions. In contrast to hadrons, they directly probe the 
entire space-time evolution of the expanding system, escaping freely without final-state 
interactions. At low masses M<1 GeV (LMR), thermal dilepton production is mediated 
by the broad vector meson p (770) in the hadronic phase. Due to its strong coupling to 
the TTTT channel and the short life time of only 1.3 fm/c, "in-medium" modifications of 
its mass and width close to the QCD phase boundary have since long been considered 
as the prime signature for chiral symmetry restoration [U El [3] . At intermediate masses 
M>1 GeV (IMR), it has been controversial up to today whether thermal dileptons are 
dominantly produced in the earlier partonic or in the hadronic phase, based here on 
hadronic processes other than vrvr annihilation. Originally, thermal emission from the 
early phase was considered as a prime probe of deconfinement [U [5] . 

Experimentally, it took more than a decade to master the challenges of very rare 
signals and enormous combinatorial backgrounds. The first clear signs of an excess of 
dileptons above the known decay sources at SPS energies were obtained by CERES [6] for 
M<1 GeV, NA38/NA50 (T] for M>1 GeV and by HELIOS-3 [8] for both mass regions 
(see |[9| for a short recent review including the preceding pp era and the theoretical 
milestones). The final status reached by CERES jTU] and NA50 [7] is illustrated in 
Fig. [H The sole existence of an excess gave a strong boost to theory, with hundreds 
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Figure 1. Excess dileptons seen in previous SPS experiments by CERES [10] (LMR, 
left) and NA50 [7| (IMR, right). For the former, there is now some sensitivity (around 
0.9 GeV) to specific theoretical predictions. 

of publications. In the LMR region, vrvr annihilation with regeneration and strong 
in-medium modifications of the intermediate p during the fireball expansion emerged 
as the dominant source. However, the data quality in terms of statistics and mass 
resolution remained largely insufficient for a precise assessment for the in-medium 
spectral properties of the p. In the IMR region, thermal sources or enhanced charm 
production could account for the excess equally well, but that ambiguity could not be 
resolved, nor could the nature of the thermal sources be clarified. 

A big step forward in technology, leading to completely new standards of the data 
quality in this field, has recently been achieved by NA60, a third-generation experiment 
built specifically to follow up the open issues addressed above [H]. Initial results 
on mass and transverse momentum spectra of the excess dimuons have already been 
published [121 113] ■ This paper shortly reviews these results, but also reports on further 
aspects associated with the centrality dependence, polarization, acceptance-corrected 
mass spectra and absolute normalization. 

Fig. [2] (left) shows the centrality-integrated net dimuon mass spectrum for 158A 
GeV In-In collisions in the LMR region. The narrow vector mesons u and are 
completely resolved; the mass resolution at the u) is 20 MeV. The peripheral data can 
be completely described by the electromagnetic decays of neutral mesons [121 [H] ■ This 
is not true for the more central data as plotted in Fig. [21 due to the existence of a strong 
excess of pairs. The high data quality of NA60 allows to isolate this excess with a priori 
unknown characteristics without any fits: the cocktail of the decay sources is subtracted 
from the total data using local criteria, which are solely based on the mass distribution 
itself. The p is not subtracted. The excess resulting from this difference formation is 
illustrated in the same figure (see |l2l [131 E] ^oi details and error discussion). 

The common features of the excess mass spectra can be recognized in Fig. [2] (right). 
A peaked structure is always seen, residing on a broad continuum with a yield strongly 
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Figure 2. Background-subtracted mass spectrum before (dots) and after subtraction 
of the known decay sources (triangles). Right: Excess dimuons compared to theoretical 
predictions, keeping the original normalization in absolute terms |15] . 



increasing with centrality (see below), but remaining essentially centered around the 
nominal p pole [l4j . Without any acceptance correction and pt selection, the data can 
directly be interpreted as the space-time averaged spectral function of the p, due to a 
fortuitous cancellation of the mass and px dependence of the acceptance filtering and 
the phase space factors associated with thermal dilepton emission [H]. The two main 
theoretical scenarios for the in-medium spectral properties of the p, broadening [2] and 
dropping mass [3], are shown for comparison. Both have been evaluated for the same 
fireball evolution [15], and the original normalization is kept (in contrast to previous 
versions of the figure [121 [II])- Clearly, the broadening scenario gets close, while the 
dropping mass scenario in the version which described the CERES data reasonably 
well [2l El E] fails for the much more precise NA60 data. A strong reduction of in- 
medium VMD as proposed by the vector manifestation of chiral symmetry [16j would 
make hadron spectral functions in hot and dense matter altogether unobservable, but 
central aspects of this scenario are totally unclear, and quantitative predictions which 
could be confronted with data have not become available up to today. 

A detailed view of the shape of the excess mass spectra is obtained by using a side 
window method fl^ to determine separately the yields of the peak and the underlying 
continuum. The left panel of Fig. [3] shows the centrality dependence of these variables: 
peak, underlying continuum and total yield in the mass interval 0.4<M<1.0 GeV, all 
normalized to the cocktail p. The continuum and the total show a very strong increase, 
starting already in the peripheral region. In the right panel of Fig. [3l the same data 
of the 271 region is plotted on a double logarithmic scale, but here the excess in a mass 
window above M> IGeV is also contained. This increases even steeper than the total 
low-mass yield, as is clearly borne out by the rising ratio of the two. The rise is about 
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Figure 3. Left: Excess yield ratios for peak, continuum and total vs. centrality for 
the mass window 0.4<M<1 GeV. Right: same as left on a double log scale, including 
here also the total in the mass window 1.0<M<1.4 GeV. 



linear implying, in view of the normalization to the cocktail p, that the absolute yield 
would be quadratic in Nch- 

The central NA60 results in the IMR region [IT] are shown in Fig. HI The use of 
the S'z-vertex tracker allows to measure the offset between the muon tracks and the 
main interaction vertex and thereby to disentangle prompt and offset dimuons from D 
decays. The left panel of Fig. H] shows the offset distribution to be perfectly consistent 
with no charm enhancement, expressed by a fraction of 1.0 ± 0.1 of the canonical level. 
The observed excess is really prompt, with an enhancement over Drell-Yan by a factor 
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Figure 4. Left: Fit of the weighted offset distribution in the IMR region with the 
contributions from charm and prompt decays. Right: Acceptance-corrected mass 
spectra of Drell-Yan, open charm and the excess (triangles). 



of 2.4. The excess can now be isolated in the same way as was done in the LMR region, 
subtracting the measured known sources, here DY and open charm, from the total data. 
The right panel of Fig. H] shows the decomposition of the total into DY, open charm 
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Figure 5. Polar angular distributions for the excess and the lu in the Collins Soper 
frame. The data are fit with the expression da / dcosOcs = 1 + acos^Bcs- 



and the prompt excess. The mass spectrum of the excess is quite similar to the shape 
of open charm and much steeper than DY; this explains of course why NA50 could 
describe the excess as enhanced open charm. The transverse momentum spectra are 
also much steeper than DY. The fit temperatures of the rriT spectra associated with 3 
mass windows are indicated on the bottom of the figure. 

The remainder of this paper is concerned with excess data fully corrected for 
acceptance and pair efficiencies [131 [IB]- In principle, the correction requires a 4- 
dimensional grid in the space of M-pT-y-cosOcs (where 9cs is the polar angle of the 
muons in the Collins Soper frame). To avoid large statistical errors in low-acceptance 
bins, it is performed instead in 2-dimensional M-p^ space, using the measured y and cos9 
distributions as an input. The latter are, in turn, obtained with acceptance corrections 
determined in an iterative way from MC simulations matched to the data in M and pt- 
The ^/-distribution is found to have the same rapidity width as dNch/drj, ay ~ 1.5 |18j . 
The cosdcs distributions for two mass windows of the excess and the uj are shown in 
Fig. [5l Within errors, they are found to be uniform, implying the polarization of the 
excess dimuons to be zero, in contrast to DY and consistent with the expectations for 
thermal radiation from a randomized system. 

The two major variables characterizing dileptons are M and pt, and the existence 
of two rather than one variable as in case of real photons leads to much richer 
information. Beyond (minor) contributions from the spectral function, p^ encodes the 
key properties of the expanding fireball, temperature and transverse (radial) fiow. In 
contrast to hadrons, however, which receive the full asymptotic fiow at the moment of 
decoupling, dileptons are continuously emitted during the evolution, sensing the space- 
time development of temperature and fiow. This makes the dilepton pj- spectra sensitive 
to the emission region, providing a powerful diagnostic tool [HI [19]. Fig. [6] (left) displays 
the centrality-integrated invariant tut spectra, where ttit = {p'^ + M'^Y^'^, for four mass 
windows; the is included for comparison. The ordinate is normalized to dNch/dr) in 
absolute terms, using the same procedure as described in detail for the cf) [20] and relating 
Npart — dNch/drj at ri=2.9 as measured to within 10% by the Si pixel telescope. Apart 
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Figure 6. Acceptance-corrected transverse mass spectra of the excess dimuons for 4 
mass windows and the (f) 13J (left), and a decomposition into peak and continuum for 
the yO-like window (right, see text). The normahzation in absolute terms is independent 
of rapidity over the region measured. For error discussion see |13j . 

from a peculiar rise at low mT (<0.2 GeV) for the excess spectra (not the 0) which 
only disappears for very peripheral collisions [9lll3], all spectra are pure exponentials, 
but with a mass-dependent slope. Fig. [6] (right) shows a more detailed view into the 
p-like mass window, using the same side-window method as described in connection 
with Fig. [3] to determine the pr spectra separately for the p peak and the underlying 
continuum. All spectra are purely exponential up to the cut-off at pr=3 GeV, without 
any signs of an upward bend characteristic for the onset of hard processes. Their slopes 
are, however, quite different (see below). 

The inverse slope parameters T^s extracted from exponential fits to the rriT spectra 
are plotted in Fig. [7] (left) vs. dimuon mass [13], unifying the data from the LMR and 
IMR regions. The hadron data for rj, uj and (f) obtained as a by-product of the cocktail 
subtraction are also included, as is the single value for the p-peak from Fig. [6] (right). 
Interpreting the latter as the freeze-out p without in-medium effects, all four hadron 
values together with preliminary vr" data from NA60 can be subjected to a simple 
blast-wave analysis. The results, plotted in a plane of freeze-out temperature Ty^ and 
average expansion velocity (/St), are shown in Fig. [7] (right). Dilepton and hadron data 
together suggest the following consistent interpretation. Maximal flow is reached by the 
p, due to its maximal coupling to pions, while all other hadrons freeze out earlier. The 
Teff values of the dilepton excess rise nearly linearly with mass up to the p-pole position, 
but stay always well below the p line, exactly what would be expected for radial flow of 
an in-medium hadron-like source (here tt+tt^^p) decaying continuously into dileptons. 

Beyond the p, however, the Tes values of the excess show a sudden decline by about 
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Figure 7. Left: Inverse slope parameter Toff vs. dimuon mass for the combined 
LMR/IMR regions of the excess in comparison to hadrons [13]. Open charm is 
subtracted throughout. For error discussion see [13] ■ Right: Blast wave results based 
on the Teff values for 7r~, rj, p, uj and (j) (see text). 

50 MeV. Extrapolating the lower-mass trend to beyond the p, such a fast transition to 
a seeming low-flow situation is extremely hard to reconcile with emission sources which 
continue to be of dominantly hadronic origin in this region. A more natural explanation 
would then be a transition to a dominantly early partonic source with processes like 
qq yU+yU^ for which flow has not yet built up [19]. While still controversial [21], this 
may well represent the first direct evidence for thermal radiation of partonic origin, 
overcoming parton-hadron duality for the yield description in the mass domain. 

The acceptance- and efficiency-corrected data can also be projected on the mass 
axis. Fig. [8] shows a set of mass spectra for some selected slices in px to illustrate 
the evolution from low to high pt- Recent theoretical results from the three major 
groups working in the field are included for comparison [191 1211 [22|. At very low pt, 
a strong rise towards low masses is seen in the data, refiecting the Boltzmann factor, 
i.e. the Plank-like radiation associated with a very broad, nearly fiat spectral function. 
Only the Hees/Rapp scenario [21j is able to describe this part quantitatively, due to 
their particularly large contribution from baryonic interactions to the low-mass tail of 
the p spectral function. At higher pt, the infiuence of radial fiow increasingly changes 
the spectral shapes, and at very high px, all spectra appear p-like. Here, only the 
Renk/Ruppert results [I9] seem to contain sufficient fiow to describe the data. 

This paper contains the most comprehensive data set on excess dileptons above the 
known sources which has so far become available through NA60. All observations can 
consistently be interpreted in terms of thermal radiation from the fireball in the whole 
mass region up to the J/ip. The superior data quality and the resulting clarity of the 
physics messages has yet to be matched by any other dilepton experiment in the field. 
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Figure 8. Acceptance-corrected mass spectra of the excess dimuons in selected slices 
of pt- Absolute normalization as in Fig. [6] 
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